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The disintegration of quasiparticle in a quantum dot due 
to the electron interaction is considered. It was predicted 
recently that above the energy e* = /S.{g / \ngY^^ each one 
particle peak in the spectrum is split into many components 
(A and g are the one particle level spacing and conductance) . 
We show that the observed value of e* should depend on the 
experimental resolution 5e. In the broad region of variation of 
5e the In g should be replaced by ln(A/(/5e). We also give the 
arguments against the delocalization transition in the Fock 
space. Most likely the number of satellite peaks grows con- 
tinuously with energy, being ~ 1 at e ~ £*, but remains finite 
at e > e*. The predicted logarithmic distribution of inter- 
peak spacings may be used for experimental confirmation of 
the below-Golden-Rule decay 

PACS numbers: 72.15.Lh, 72.15.Rn, 73.23.-b 



Decay of single-electron excitations in quantum dots 
became now the subject of intensive experimental 
and theoretical |^|-^ investigations. For a closed quan- 
tum dot instead of real decay of a quasiparticle one 
should consider the disintegration of one (5-peak in the 
single-particle spectral density p[e) into a relatively dense 
bunch of peaks. Each component of this bunch repre- 
sents the one particle contribution into complicated ex- 
act eigenstate. Quasiparticle life-time in large Fermi sys- 
tem is usually associated with decay into two-particle- 
one-hole configuration. The corresponding width may 
be found using the usual Golden Rule . The energy e' 
at which this width becomes of the same order of magni- 
tude with three-particle level spasing gives us the natural 
threshold for decay of a quasiparticle. However, it was 
shown by the authors of Ref. [|j that due to the effec- 
tive interaction with five-particle, seven-particle and so 
on excitations (they call all states consisting of n + 1 par- 
ticles and n holes 2n -I- 1-st generation) the actual thresh- 
old for particle disintegration is much lower. The more 
detailed investigation of statistics of states constituting 
the single-particle excitation below e' is the subject of 
this letter. The statistical approach to finite interacting 
Fermi systems has a long history [[7[-pl|. However, the 
main attention was paid to the investigation of the fully 
developed chaos. Here we are interested in the very be- 
ginning of the chaotic behaviour, then quasiparticle may 
be coupled with only a few many particle states and spac- 
ing between peaks which we consider is much larger than 
the level spacing. 

The convenient quantity, which describes the splitting 
of a noninteracting quasiparticle peak into many peaks. 



is the participation ratio (PR) P — ctf ■ Here of is the 
relative strength of individual peak in p{e) and the sum 
over bunch of peaks corresponding to one particle exci- 
tation is ^ af = 1 . Physically the PR is the inverted ef- 
fective number of exact many particle eigenstates consti- 
tuting the quasiparticle excitation. From technical point 
of view, the authors of Ref. Q have summed up starting 
from the small excitation energies the series of special 
perturbative contributions leading to quasiparticle dis- 
integration and then estimated, at which energy e* this 
series blows up. In terms of PR this procedure gives 
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where A is the averaged single-particle level spacing, 
5 ^ 1 is the dimensionless conductance, e ^ A is the 
energy of our quasiparticle, and Pn are some numerical 
coefficients. Each new term in the sum in Eq. (^ cor- 
responds to taking into account more and more compli- 
cated admixtures to quasiparticle. The first term {n = 0) 
describes the mixing with 2-particles and 1-hole, second 
to 3-particles and 2-holes, and so on. 

At energies close to e* ~ ^{g/ Iny)^/^ all terms of the 
series in Eq. (|]) become of the same order of magnitude. 
This means that at e > e* the PR can not be close to 
1. However, the concrete way of the quasiparticle disin- 
tegration with the growth of e depends on the behaviour 
of the coefficients p„. In general, the three possibilities 
for the asymptotics of p„ are 



a). Pn 



■,b). Pn ^ aJ^n' ■,c).pn^l/n\ 
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which corresponds to zero, finite and infinite radius of 
convergence of the series in Eq. (^ respectively. One 
should naturally expect very different features of the re- 
sumed result in these three cases. Mapping the prob- 
lem of quasiparticle lifetime onto that of particle hopping 
on the Cay ley tree led the authors of Refs. to the 
asymptotics (^). However, as we will show below, the 
actual asymptotics is close to the Eq. (j^c). 

In general, in order to observe experimentally the split- 
ting described by Eq. (|l|) one should be able to re- 
solve all many particle eigenstates, which means that 
the experimental errors should be exponentially small 
5e ~ Aexp(-27r^e/6A) |l|]. Therefore, in this letter 
we would like to find how the mechanism considered in 
will manifest itself for more realistic de. First of all, 
even in order to see the decay of quasiparticle into three- 
particle configurations one should have sufficiently good 



resolution 5e ~ A^/e^ (any few peaks falling into the seg- 
ment Se are seen as one of joint strength J^Ss '^i)- 
most interesting is the result for A^/e"* <C ^ A'^/e^. 
Physically this means that accuracy is much better than 
needed to resolve the three-particle levels, but not enough 
to see the five-particle ones. In this case 

^ = 1 - -^bie) , bis) = ^ 6„ , (3) 

where Sc = A 1/17/ ln(A/gfc). The transition from pure 
single-particle to split spectrum now takes place at e ~ 
Ec- In particular if Se ~ A^/e^ one has £c ~ e' = Ay/g 
in accordance with the Golden Rule prediction |^,^. At 
Se ^ /e^ the expansion formally is restored, but 
the coefficients of this new series p* are much smaller 
than those of the Eq. (|l|). For better accuracy Se <C 
A^/e^ the coefficients p* become a functions of the reso- 
lution p* — p*^{Se). Only at extremely small Se one has 
A"+Ve")«p„. 
As it was shown in Refs. [|[^ the values of the matrix 
elements (MEs) of two-particle interaction are Gaussian 
distributed with the variance: 

T^=AVg' . (4) 

Here the numerical factors ^ 1 (see e.g. Ref. [||) are 
included into the definition of 17 ^ 1 The estimate 
(^ was done for the diffusive quantum dot. However, 
our approach may be valid for the ballistic dot also. The 
only necessary condition is that the MEs of interaction 
should be random with the amplitude \V\ ^ A. 

Gonsider first the mixing of particle with three-particle 
states (two particles and one hole). The density [dn/de) 
of this statesNo spin of quasiparticle. is 

i^3=£V4A3 ■ (5) 

Here one factor 1/2 comes from the integration over the 
three energies at fixed epi + ep2 + en = e and another is 
added due to the Fermi statistics of two produced par- 
ticles. We are interested in energies e ^ A^. There- 
fore \V\v3, ^ 1, which means that the majority of one- 
particle states are almost nonperturbed. In this case the 
main contribution to the PR comes from the very small 
part of levels, for which the energy difference between 
one- and three-particles excitations e*^^-' — e'^^^ occasion- 
ally turns out to be of the same order of magnitude with 
the ME V . The relative fraction of such states is small 
~ but their PR differs by 100% from P = 1. 

Therefore, the averaged contribution of such events to 
P is SP3 ~ \V\v^ ~ e^/^A^. The accurate calculation 
allows to find also the numerical factor 

F3 = 1 - 27rjTn'K3 = 1 - yJVj^e^lgA^ . (6) 

Here due to the Eq. (|) M = ^J^^/g- 



Mixing of quasiparticle with higher generations (five- 
particle, seven-particle etc.) may be formally taken into 
account in the same way 

5P2n+i = -2n\V^y;^^\y2n+i , (7) 

where i'2n+i ^ j ^n+i ^ rp-j^g only difference from 
(^ is that now vj^^'''^^ is the effective ME connecting 
the first and 2n -I- 1-th generations via the n-th order 
of the usual perturbation theory. The naive estimate of 
this effective interaction gives 14// ~ (A/(7)"(l/A)"~-'^, 
which lead to SP2n+\ ^ —{e^ j gt^Y' . The main ad- 
vantage of Ref. jij was in fact the observation that the 
high order corrections to P have additional enhancement 
~ (In g)"^^ compared to the naive estimate. In order to 
demonstrate the origin of this large logarithms, consider 
the effective ME connecting generations 1 and 5 

iV^ ^^121^23 , 2A2 de 2A 
V/f = = — 2] TT — jlnff-(8 

Here we have left in the sum over £2 only one level closest 
to £1 (contribution of the other levels is ~ 1/ In g smaller) 
and then averaged over its position. Therefore, the upper 
bound of the integral is |£| < A. More interesting is 
the origin of the lower bound. The use of the effective 
interaction requires |Vi2|, IV23I < |£i - £2! ~ |£3 - £2]- 
Otherwise one should consider the strong mixing of three 
almost degenerate states |1), |2), |3) (accurate taking into 
account such a 3-level interaction leads also to ~ 1 / In g 
corrections) . That is why the lower bound in the integral 
in Eq. (1) is |£| > A/g. 

The ME (1) together with vr^ ^ e^ / A^ allows one to 
reproduce the first nontrivial term of the expansion (|^). 
Moreover, even if one does not take into account the mix- 
ing with higher generations, the correction described by 
(|) could blow up the PR (|) at £ ~ Ag^''^{hig)-''''^ . 
However, one more important feature of p{e) may be 
demonstrated by the Eq. (|8|). The logarithmic diver- 
gence of the integral in (|^) shows that the MEs with 
very different denominators are equally important for the 
PR. Suppose we are not able to resolve peaks which are 
closer than some Se <^ A/g. As we mention before, one 
can see two peaks of comparable amplitude at distance 
Se only if Se ~ T4// ~ £1 — £3. This means that the 
upper bound in the integral in Eq. (||) should be chosen 
via Veff ~ A^/(7^£i2 > Se and thus 

._^|j,^l^^^^ln(A) , ,„ 

in accordance with (||). In order to illustrate this re- 
sult we have shown on the Fig. 1 the density of cou- 
ples of peaks as a function of the logarithm of spacing A 
between them dn/dln{X) (at e slightly below £*). The 
mixing with third generation leads to the narrow (width 



2 



^ 1) peak at A ~ A'^/e^ ~ ^/ff- The contribution from 
fifth generation at A '--^ A/g is in e^/gA^ times weaker, 
but such events are uniformly distributed over the wide 
region IniA/g"^) < ln(A) < ln(A/g). 

GeneraHzation of (p|) for an arbitrary generation gives 
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The upper limit for all integrals here is the same as in 
£i < A. The small values of Ei are restricted due to 
the requirement that none of the intermediate states in 
Veff could be mixed strongly with initial or final state. 
One may found the lower and upper bounds for |Ve//| 
by considering the simplified version of the logarithmic 
inequalities in (|l0|): \n{gei/A) > for any i for lower 
bound and > for upper bound. For 

large n such calculation gives 



1<A< 



(11) 



Thus at least the integral (|T^) could not contain any n! 
lH. Eqs. (0jlOl) together lead to the Eq. ([|). 

For finite accuracy one should take into account only 
the MEs exceeding the experimental error Veff > Se, 
which is equivalent to the additional restriction on the 
domain of integration 



^ln(5e,/A) <ln(A/.g&) 
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If in addition \n{A/gde) ^ the integration in Eq. 

( p^ ) may be performed explicitly 
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In terms of log-distribution of level spacings dn/d\ii{X) 
shown on the Fig. 1 the contribution of generation 2n + 1 
leads to correction ~ [ln(A/(7(5e)]"^^. 

Consider now the physical consequences for the spec- 
trum of the different variants of asymptotic behaviour of 
the coefficients p„ (6„) shown in Eq. (^: 

a) . In fact, there is no real danger in divergence of the 
asymptotic series. One should simply break the summa- 
tion at the smallest term (with the number Uc ~ \J e* je 
or ric ~ \/ Sc/s) . The same smallest term gives the order 
of magnitude estimate of the rest(nonperturbative) part 
of the sum. p{e),b{e) become completely nonperturba- 
tive at e > £*, Sc- However, there is now indication that 
it should be p{e > e*) ^ 1. One may even see no consid- 
erable disintegration of quasiparticle peaks at £ < A^/g. 

b) . The series in £^ has finite radius of convergence 



R = £*^/a {si /a) and the 7 is responsible for the 
kind of singularity of the resummed result at £^ = i? 
(both a, 7 ^ 1). Close to this point all terms of the se- 
ries become equally important. It is natural to consider 
such a behaviour as an indication of the localization- 
delocalization transition in the Fock space |^,|| . 
c). The series is absolutely convergent. We consider this 
as the indication of absence of delocalization transition. 

The estimates of ME (||,|l^,|l^ were done for a given 
tree-type Feynman diagram connecting given initial and 
final states. Now we have to estimate the number of such 
diagrams. First of all, the density of final states: 
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Here (2n)! appears after the integration over energies of 
final particles(holes), {n+l)\ and n\ account for the 71 + 1 
identical particles and n holes. The number of diagrams 
for fixed final state is easy to estimate for the Schrodinger 
perturbation theory. The examples of diagrams for the 
screened Coulomb interaction V{x — y) ~ 5{x — y) are 
shown in Fig. 2 Each individual ME oiV{x — y) 

corresponds to decay of one particle into two particles 
and one hole, ore one hole into two holes and one particle. 
In order to find the number of diagrams it is convenient 
to start from the final state. At first step there are (n -|- 
l)n'^ /2 ways to join two of {n + 1) particles and one of 
n holes into one particle and (n -I- 1)71(7?, — l)/2 ways to 
join one particle and two holes into one hole. Then the 
same procedure may be repeated with n particles and 
77 — 1 holes. The number of diagrams connecting the 
same initial and final states found in this way is 



2-"n!(r7+ l)!(2r7- 1)!! 



(15) 



The doubling of single-particle peaks is based on very 
rare events of almost coincidence of the small ME and 
small energy difference. This means that the probability 
to find two equally large MEs is small and one should 
simply multiply the correction (^ by the number of sta- 
tistically independent diagrams. However, not all of the 
diagrams ( p^ are statistically independent. First of all, 
we have not taken into account the Fermi statistics of 
particles in the intermediate states. This means, that 
some of the diagrams should cancel each other. Second, 
we have estimated the number of diagrams of Schrodinger 
perturbation theory. If one goes to the Feynman technic, 
many of the diagrams having the same MEs and different 
energy denominators will be joined into one. For exam- 
ple, for two diagrams of Fig. 2 one has 



l/£a£b + l/£a£c = l/eh£c 



(16) 



because Sb + £c = £a (for almost degenerate initial and 
final states). Here £a,b,c are the energy denominators for 
corresponding cross section on the figure. Therefore, the 
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Eq. (15) gives only the upper bound of the number of 
independent diagrams. Combining together ( |l^ , p^ and 
the estimate of \Veff\ one finds 



Pn < {const) 



„n!(n + l)!(27i- 1)!! (const)" 



(2n)!(n+ 
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We see that combinatorics of the diagrams (jl^) could not 
compensate the decrease of phase space and the asymp- 
totics of Pn (as well as 6„) is described by the Eq. (^). 
Slightly above e — e*,£c due to the mixing with finite 
number (^ \/ln g) high order (with n ~ In 5) genera- 
tions the PR becomes sufficiently smaller than one. This 
finite number of connected generations constitutes the 
main difference of our result from what happens on the 
Cayley tree [Hjs), where even the first splitting of the 
quasiparticle peak into two proceeds through the inter- 
action with all generations. For higher energies our per- 
turbative approach formally is not valid. We are able to 
consider rigorously only the first splitting of quasiparti- 
cle peak into two. In order to go further one should be 
able to diagonalize exactly the three-levels almost degen- 
erate events, then the four-levels and so on. Mathemat- 
ically, this means that one has to sum up the series of 
~ l/\n{g) corrections to P. Nevertheless, it is natural 
to suppose, that further disintegration of quasiparticle 
also proceeds through the interaction with finite number 
of generations. If so, the number of peaks constituting 
one excitation most likely will grow smoothly with energy 
(crossover instead of phase transition). The delocaliza- 
tion in the Pock space will not take place in this scenario 
(although, it may be difficult to find the experimental 
evidence of presence or absence of such delocalization). 

Even more informative than P is the distribution of 
spacings inside the quasiparticle bunch. The distribu- 
tion of spacings for first decay into two peaks (two dis- 
tinct bunches) has complicated hierarchical structure. 
The natural variable to describe this distribution is In A 
(Pig. 1). In particular this means that the disintegration 
threshold Ec should depend on the experimental accu- 
racy. It is natural to expect that this log-distribution 
of spacings will survive after further disintegration into 
three and more peaks. Moreover, both new delocaliza- 
tion thresholds e* and Ec differ only by the square root 
of the logarithm from the Golden Rule prediction, which 
makes them quite difficult to be observed in the direct ex- 
periment. However the wide logarithmic distribution of 
spacings within the single particle bunch of peaks (like 
that on the Pig. 1) may be easily distinguished from 
e.g. Poisson or Wigner-Dyson distribution. Thus we 
may conclude that the investigation of spacings distri- 
bution in the single particle spectral density should open 
the easiest way to observe the below-Golden-Rule decay 
of quasiparticles in quantum dot predicted in Ref. 
Also, the further investigation of quasiparticle decay may 
be performed numerically. 
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FIG. 1. Distribution of spacings A for first doubling of the 
peaks as a function of In A. The mixing with generations 
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FIG. 2. The examples of diagrams. Energy denominators 
are associated with transverse sections (dashed lines). 



